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ABSTRACT 

We present UBVRI CCD photometry in the region of the open clusters NGC 637 
and NGC 957. The radii are found to be 4.'2 and 4.'3. Their reddenings E(B - V) 
are 0.64 ± 0.05 mag and 0.71 ± 0.05 mag and their distances, from main sequence 
fitting are 2.5 ±0.2 kpc and 2.2 ±0.2 kpc. Comparison with Z = 0.02 isochrones leads 
to an age of 10 ± 5 Myr for both clusters. Combining our photometry with 2MASS 
JHK shows the reddening law in these directions to be normal. Mass function slopes 
of x = 1.65 ± 0.20 and 1.31 ± 0.50 are derived for the clusters, both of which are 
found to be dynamically relaxed. Spectral and photometric characteristics of three Be 
stars, two in NGC 957 and one (newly discovered) in NGC 637 indicate them to be of 
Classical Be type." 

Key words: Key words: stars: emission-line, Be HertzsprungRussell (HR) diagram 
stars: luminosity function, mass function ISM: dust, extinction Galaxy: open cluster 
and associations: individual: NGC 637 and NGC 957. 



1 INTRODUCTION 

Open clusters are excellent targets to understand issues re- 
lated to Galactic structure, chemical composition, stellar 
population, dynamical evolution and star formation pro- 
cesses in the Galaxy. The importance of photometric stud- 
ies of star clusters lie in their colour-magnitude diagrams 
(CMDs) which allow us to estimate cluster's fundamental 
parameters such as age, distance and reddening. In this 
study, we consider two young open cluster, NGC 637 and 
NGC 957, situated in perseus arm of Milky Way to estimate 
their fundamental parameters and mass function. 

NGC 637 = C0139+637 (q 20 oo = 01 h 43 m 04 s , 
<5 2 ooo=+64°02'24"; I = 128? 55, b=l°73 ) is classified as a 
Trumpler class I 2m by iRuprechtl (Il966t ). Using RGU pho- 
tographic data, iGrubissich] (Il975l ) found the cluster to 
be located in the Perseus arm at a distance of 2.4 kpc. 
iHuestamendia et al.l ()l99ll ) used the UBV photoelectic pho- 
tometry and derived a colour excess of E(B — V) = 0. 6 6 mag 
and a distance of 2.5 kpc. Recently, IPhelps fc Janesl (|l994l ) 
studied this cluster using UBV CCD photometry and es- 
timated a colour excess of E(B — V) = 0.49 ± 0.03 mag. 
Further, using CMD and assuming a total-to-selective ab- 
sorption (Rv) of 3.1, they derived a distance of 2.75 kpc. 

* E-mail: rkant@aries.emet. in 



NGC 957 = C0230+573 (q 20 oo = 02 h 33 m 21 s , 
5 2 ooo =57° 33' 36"; J=136?28, 6=-2?65) is locat ed in the 
Perseu s spiral arm of our Galaxy and according to lRuprechtJ 
dl966h. it is of type III 2 p. Th is cluster was studied by 
ICimenez fc Garcia-pelavol (jl980l l using RGU photographic 
photometry and they found that this is an young star cluster 
situated at a distance of 1.8 kpc. T his cluster was further in- 
vestigated bv lGerasimenkol (|l99ll ) using Johnson UBV pho- 
toelectric and photographic photometry. He derived redden- 
ing of E(B — V) = 0.90 mag, distance of 2.1 kpc and age of 
about 4 Myr. 

In the present study, we provide deep UBVRI CCD 
photometry of NGC 637 and NGC 957, and re- visit their ba- 
sic parameters and estimate their luminosity and mass func- 
tions. Both the clusters are poorly studied in the literature. 
For NGC 957, CCD photometric observations is presented 
for the first time alongwith spectroscopic observations of a 
few brightest members of both the clusters. 

The paper is organized as follows. We present the obser- 
vations and data reductions in Sj2] while in £j3j we describe 
the data analysis. Sj4]and Sj5]describe about the mass function 
and mass segregation. Finally, $6]contains the conclusions of 
our study. 
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Figure 1. Finding chart of the stars in the field of NGC 637 and NGC 957. The (X, Y) coordinates are in pixel units corresponding to 
0".72 on the sky. Direction is indicated in the map. Filled circles of different sizes represent brightness of the stars. Smallest size denotes 
stars of F~20 mag. Open circle represent the cluster size. SI star in the region of NGC 637 is emission type star. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometric observations 

Johnson UBV and Cousins RI CCD photometric observa- 
tions were obtained using 104-cm Sampurnanand telescope 
at Nainital on 13 November 2004 for NGC 637 and on 28 
September 2006 for NGC 957. Data were recorded using 
a liquid nitrogen cooled 2kx2k CCD of 24/xm square size 
pixel, which results in a scale of / .'36 per pixel and a square 
field of 12!7 in the sky on a side. The CCD has a gain of 



10 electrons per analog-to-digital unit and a read out noise 
of 5.3 electrons. To improve the signal-to-noise ratio, all the 
observations were taken in 2 x 2 binning mode. Log of obser- 
vations are given in Table 1. In V and I filters, we also ob- 
served field regions situated about 15'north from the cluster 
center of both NGC 637 and NGC 957. Identification maps 
of cluster and field regions are shown in Fig. 1. 

A number of bias and flat-field frames were taken dur- 
ing the observations. Flat-field exposures were taken of the 
twilight sky in each filter. Bias and flat field corrections were 
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Table 1. Log of observations, with dates and exposure times 
for each passband. N denotes the number of stars measured in 
different passband. 



Band 


Exposure Time 
(in seconds) 


Date 




N 






NGC 637(Cluster 


region) 




U 


1800X2, 300x1 


2004 November 


13/14 


1057 


D 


1200x3, 300x1 


2004 November 


13/14 






120x1; 10x1 


2004 November 


13/14 


2392 


V 


900x3, 120x1; 10x1 


2004 November 


13/14 


3261 


R 


480x3, 60x2; 5x1 


2004 November 


13/14 


3219 


I 


300x3, 60x1; 5x1 


2004 November 


13/14 


3151 






NGC 637(Field region) 




V 


180x2 


2004 November 


14/15 


1286 


I 


180x2 


2004 November 


14/15 


1266 






NGC 957 






u 


1800X2, 300x1 


2006 September 


28/29 


915 


B 


1200x2, 300x2 


2006 September 


28/29 


2087 


V 


900x2, 180x2 


2006 September 


28/29 


2855 


R 


600x2, 120x2 


2006 September 


28/29 


2819 


I 


300x3, 60x2 


2006 September 


28/29 


2707 






NGC 957(Field region) 




V 


180x2 


2004 November 


14/15 


1284 


I 


180x2 


2004 November 


14/15 


1269 




14.0 16.0 18.0 20.0 22.0 

V mag 



performed using the standard tasks. Subsequent data 

reduc tion and analysis we re done using the DAOPHOT soft- 
ware (|Stetsonlll987l . ll992l V We performed the profile fitting 
photometry using quadratically varying point spread func- 
tion. For each filter, the stars have been aligned to that of 
a reference frame with largest exposure time and an aver- 
age instrumental magnitude weighted by the photometric 
error was derived. A final catalogue was created of stellar 
objects identified in at least two filters, with shape-defining 
parameter —2 < sharpness < 2 and and the goodness-of-fit 
estimator x < 5. 

For translating the instrumental magnitude to the 
standard magnitude, we observed the standard field PG 
0231+051 (|LandolrJ \l99% ) on 20 November 2004 several 
times during the night to determine the Earth's atmospheric 
extinction coefficients, color coefficients and the zero points 
of the telescope and CCD system. The standard stars used in 
the calibrations have brightness range of 12.77 < V < 16.11 
and a color range of -0.329 < (B-V) < 1.448. The transfor- 
mation equations derived using linear least square regression 
are as follows: 

u = U + 6.92 ± 0.01 - (0.06 ± 0.01) (C/ -B) + (0.54 ± 0.02) A 
b — B + 4.66 ± 0.01 - (0.02 ± 0.01)(B - V) + (0.24 ± 0.01)X 
v = V + 4.22 ±0.01 - (0.01 ± 0.01) (B-V) + (0.14 ±0.01) X 
r = #±4.13±0.01 - (0.00 ± 0.01)(V - 7?) ± (0.10 ±0.01)X 
i — I + 4.60 ± 0.01 - (0.02 ± 0.02) (7? - I) + (0.06 ± 0.01)A 

where u,b,v,r and i are the aperture instrumental magni- 
tudes and U, B, V, R and I are the standard magnitudes and 
X is the airmass. Zero-point and colour-coefficient errors are 
~ 0.01 mag. 

The internal errors, as derived from DAOPHOT, in 
magnitude and colour are plotted against V magnitude in 
Fig. [21 This figure shows that photometric error is < 0.01 



Figure 2. The photometric errors (in magnitude) corresponding 
to the brightness measurement at V . U — B, B — V, V — R and 
V — I, are plotted against the V-band brightness. Errors on the 
y-axis represent the internal errors as estimated by DAOPHOT 
routine. 



Table 2. Comparison of our photometry with others for NGC 637 
and NGC 957. The difference (A) is always in the sense present 
minus comparison data. The mean along with their standard de- 
viations in magnitude are based on N stars. Few deviated points 
are not included in the average determination. 



< av > 

Jean±<r(N) 



< A(fl - V)> 
Mean±(r(N) 



A(C7 — B) > 
Mean±cr(N) 



NGC 637 



10.0 - 11.0 
12.0 - 13.0 
13.0 - 14.0 
14.0 - 15.0 



0. 01±0. 01(2) 
0.03±0.05(3) 
0.05±0.04(8) 
0.02±0.02(7) 



at al. (1991) 



Phelps & Ja 



with 
(1994) 



10.0 
12.0 
13.0 
14.0 
15.0 
16.0 
17.0 
18.0 
19.0 



11.0 
13.0 
14.0 
15.0 
16.0 
17.0 
18.0 
19.0 
20.0 



12.0 - 13.0 
13.0 - 14.0 
14.0 - 15.0 



0.05±0.06(4) 
0.05±0.04(2) 
0.08±0. 06(12) 
0.06±0. 07(13) 
0.04±0. 04(12) 
0.03±0. 04(44) 
0.03±0. 05(61) 
0.03±0. 08(67) 
0.03±0. 07(102) 

NGC 957 



.cz & Garcia-pelayo 
(1980) 

-0.06±0.09(5) 
-0.04±0.07(6) 
-0.03±0.06(3) 



0. 01±0. 01(2) 
0. 01±0. 01(3) 
0. OliO. 02(8) 
0. OliO. 01(7) 



0.03i0.01(3) 
0.02i0. 04(14) 
0.03i0. 03(15) 
0.01i0.05(2l) 
0.02i0. 04(45) 
0. OliO. 05(65) 
0. OliO. 06(73) 
0. OliO. 09(96) 



-0.03i0.05(5) 
-0.03i0.06(6) 
-0.03i0.06(3) 



0.02i0.02(2) 
0.03i0.02(3) 
0.03i0.04(8) 
0.04i0.02(4) 



0.13i0.03(3) 
0.13i0. 06(14) 
0.14i0. 08(14) 
0.13i0. 05(15) 
0.12i0. 06(27) 
O.lOiO. 07(42) 
O.lliO. 08(27) 



0.03i0.05(5) 
0.03i0.06(6) 
0.03i0.06(3) 



1 1RAF is distributed by the National Optical Astronomical Ob- 
servatory which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Sci- 
ence Foundation 
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Table 3. Log of slitless and slit spectroscopic observations, with 
dates and exposure times. 



Object 


Band 


Exposure Time 
(s) 


Date 


NGC 637 


R 


5 


2004 January 28 




R+Gr5 


60 


2004 January 28 




R+Gr5 


600 


2004 January 28 


NGC 957(East) 


R 


5 


2004 January 29 




R+Gr5 


60 


2004 January 29 




R+Gr5 


600 


2004 January 29 


NGC 957(Wcst) 


R. 


5 


2004 January 29 




R+Gr5 


60 


2004 January 29 




R+Gr5 


600 


2004 January 29 


NGC 957(West) 


R. 


5 


2004 January 29 




R+Gr5 


60 


2004 January 29 


NGC 637(1) 


Grism 7/1671 


900 


2006 September 30 




Grism 8/1671 


900 


2006 September 30 


NGC 957(1) 


Grism 7/1671 


300 


2006 September 30 




Grism 8/1671 


300 


2006 September 30 


NGC 957(2) 


Grism 7/1671 


600 


2006 September 30 




Grism 8/1671 


600 


2006 September 30 



20 22 
V 



Figure 3. Comparison of present CCD UBV photometry with 
that available in the literature. For NGC 637, th e solid (•) circles 
represent data from [Hucstarncn dia et al.l Il 991) while the open 
circles (o) come from lPhelps fc Janes! Jl994h . For NGC 957, the 
cross es (x) represent com parison with photographic data taken 
from iGerasimenkol 1 119911) . 



mag at V ~ 19 mag while error in different colours is < 0.02 
mag. The final photometric data are available in electronic 
form at the WEBDA siteQ and it can also be obtained from 
the authors. 

We cross identified our objects with that present in the 
literature. Fig.[3]shows the plots of difference A in V , (B—V) 
and (U — B) with V magnitude. Solid and open circles rep- 
resent the difference with photoelectric and CCD data. The 
average difference alongwith their standard deviation are 
listed in Table 2. The comparison of present and photoelec- 
tric data in V magnitude and (B — V), (U — B) colours in 
NGC 637 do not show any significant difference with stel- 
lar magnitudes. Comparison of NGC 957 CCD data with 
photographic data also do not show any significant differ- 
ence. So, we conclude that the present CCD photometry is 
in good agreement with photoelectic and photographic data 
for NGC 637 and NGC 957 respectively. O n the other hand 
for NG C 637, comparison of CCD data with lPhelps fc Janes! 
l| 19941 ) show a zero-point offset of ~ 0.04 mag in V and ~ 
0.12 mag in (U — B) while there is no noticeable differ- 
ence seen in (B — V). A possible cause of the zero-point 
offs et may be due to the difference in dates of observation 
for lLandoltl (|l992T i standards and the target cluster fields. 
Since the zero-points depend on the factors, such as ambient 
temperature and Earth's atmospheric conditions, these coef- 
ficients change from night-to-night. A typical nightly varia- 
tion in zero-points is found to be ~ 0.08 mag in U and ~ 0.05 
mag in B, V, R and I passband. In addition to this ~ 0.02 
mag nightly variation is found in the extinction coefficients. 



2.2 Spectroscopic observations 

We obtained slitless spectra of stars in the cluster regions 
using Himalayan Faint Object Spectrograph and Camera 
mounted with the 2-m Himalayan Chandra Telescope at 
Hanle, Indifl Grism 5 (A ~ 5200 - 10300 A, A/AA ~ 870) 
was used as dispersing element in combination with R fil- 
ter (A ~ 7000 A, AA ~ 2200 A). The slitless spectra were 
recorded on a 2kx2k CCD with 15pim pixel size and with 
an image scale of 0'.'297 per pixel, which covered a total 
area of around 10' x 10'. A typical spectral resolution of 

6 A is achieved at H a . We list the log of observations in 
Table 3. We obtained the slit spectra of the three identi- 
fied emission line stars on 30 September 2006 using Grism 

7 (3800 — 6800A) and 167 /im slit combination in the blue 
region which gives an effective resolution of 1330. The spec- 
tra in the red region were obtained on the same night using 
Grism 8 (5800 - 8350 A) and 167 /an slit, which gives an ef- 
fective resolution of 2190. The Gaussian FWHM seeing was 
about l'.'5 (~ 5 pixel). 



3 ANALYSIS OF THE DATA 

3.1 Cluster radius and radial stellar surface 
density 

We perform star counts in concentric rings around an es- 
timated centre of the cluster, and then divide it by their 
respective areas. The cluster center is estimated iteratively 
by calculating average X and Y position of the stars within 
400 pixels from an eye estimated center, until they converged 
to a constant value. The cluster center in pixel units derived 
by this method are (515, 515) for NGC 637 and (430, 630) 
for NGC 957. The cluster center in celestial coordinate are 
(02000 = 01 h 43 m 04 a .2, 5 20 oo = 64°02'23'.'5) for NGC 637 
and (aaooo = 02 h 33 m 23 s .9, <5 2 ooo = 57°33'29'.'9 for NGC 
957. The derived center of both the clusters are close to the 
one given in WEBDA (*JTJ) . In Fig. [4] we show the radial 
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Figure 4. Surface density distribution of stars in the field of 
the cluster NGC 637 and NGC 957. Errors are determined from 
sampling statistics(=l/V r /V where N is the number of stars used 
in the density estimation at that point). Vertical and horizontal 
arrows represent the radius of the clusters and field star density 
respectively. 



except in V, (U — B) CMDs where it is only up to V ~ 
19 mag. Main-sequence (MS) of the clusters is clearly vis- 
ible upto V ~ 17.0 mag. The MS fainter than 17.0 mag 
has more scatter and field star contamination is also more 
evident for fainter stars and because of this it is hard to 
delineate the field stars from the cluster members, only on 
the basis of their closeness to the main populated area of the 
CMD as field stars at cluster distance and reddening also oc- 
cupy the same area. Proper motion or radial velocity data 
are very useful in separating cluster members from the field 
stars. But due to unavailability of this kind of data in the 
literature we use photometric criterion to separate probable 
cluster members from the field stars. We selected members 
by defining the blue and red envelope around the MS which 
is shown in the CMDs of both the clusters. A star is con- 
sidered as a non-member if it lies outside the marked region 
in the CMDs. In Table [4] we have listed the expected num- 
ber of field stars using V , (V — I) CMD of the field region. 
From this Table we can estimate the frequency distribution 
of stars in different parts of the CMD. It is also clear that 
all photometric probable members cannot be cluster mem- 
bers and non-members should be subtracted in the studies 
of cluster mass function etc. However, probable members lo- 
cated within a cluster radius can be used to determine the 
cluster parameters, as they have relatively less contamina- 
tion due to field stars and the same has been assumed in the 
subsequent analysis. 



density profile for both the clusters which is calculated us- 
ing successive annuli of 50 pixels (~ 35") around the cluster 
centre. For NGC 637, it flattens around r = 250"and begin 
to merge with the background stellar density. Therefore, we 
consider 250"as t he cluster r adius. A smooth dashed line 
represents a fit by King (1962) profile: 



fir) 



fo 



1 + (r/r c )- 



(1) 



where fo is the central density, r c is the core radius of the 
cluster and /j, is the background density. For NGC 637, these 
parameters are fo = 10.96 x 10 -3 star/arcsec 2 , r c = 39'.'6 
and f b = 6 .15 x 10~ 3 star/arcsec 2 . In case of NGC 957, 
iKind l|l962l ) profile fitting could not converge because of 
large errors in the value of f(r). Radius of this cluster is 
assumed as 260"considering the field star density on the 
basis of last two points in the radial density profile. For NGC 
637 an d NGC 957 , the radius has been listed as 100"and 
300 in iDias et al.1 (|2002l ) catalog. Present estimate of the 
radius is larger for the cluster NGC 637 while it is smaller 
for NGC 957. 



3.2 Colour-magnitude diagrams 

The V, (U - B) ■ V, (B - V) and V,(V-I) CMDs of NGC 637 
and NGC 957 alongwith the V, (V - I) CMD of the corre- 
sponding field regions are shown in Fig. [5] To get the clear 
sequence of the cluster we plot only the stars within the 
cluster radius. It is evident that the morphology of CMDs 
for NGC 637 and NGC 957 is typical of a young age open 
star clusters. The CMDs extends down to V ~ 21.0 mag 



3.3 Colour-colour diagram 

To estimate the interstellar extinction towards the clusters 
we plot (U — B) versus (B — V) diagrams using probable clus- 
ter mem bers. The intrin s ic zer o-age main-sequence (ZAMS) 
given bv lSchmidt-Kalerl (|l982r i is fitted by the dotted curve 
assuming the slope of reddening E(U — B)/E(B — V) as 
0.72. The ZAMS fitted to the MS stars of spectral type ear- 
lier than AO provides a mean value of E(B — V) = 0.64±0.05 
mag for NGC 637 and E(B-V) = 0.71 ±0.05 for NGC 957. 
Our derived values of reddening agree fairly well with the 
values estimated by others (52). 



3.4 Interstellar extinction in near-IR 

2MASS JHK data is used to study the interstellar ex- 
tinction in combination with the optical data. The K s 
magnitude are con verted into K magnitude following 
iPersson et ail l|l99Sl ) . The (J — K) versus (V - K) diagrams 
for b oth the cluster s are s hown in Fig. 7. The ZAMS taken 
from ISchaller et all (|l992T l for Z = 0.02 is shown by solid 
line. The fit of ZAMS provides E(J - K) = 0.23 ± 0.20 
mag and E(V - K) = 1.70 ± 0.20 mag for NGC 637 and 
E(J - K) = 0.29 ± 0.20 mag and E(V - K ) = 2.00 ± 0.20 
mag for NGC 957. The ratios ^Iv-h] ~ 0.14 ± 0.30 and 
0.15 ± 0.20 for NGC 637 and NGC 957 respectively are in 
good agr eement with the norm al interstellar extinction value 
of 0.19 l|Cardelli et al.lll989h . However, scattering is more 
due to large error in JHK data. 
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Figure 5. The V, (U - B); V, (B - V) and V, (V - I) CMDs for the cluster NGC 637 and NGC 957 using stars within the cluster 
radius and V, (V — I) CMDs for the corresponding field regions. Black triangles show the emission stars observed by us. Solid lines 
represent the blue and red envelope of the cluster MS. 



3.5 Extinction Law 

To investigate interstellar extinction law towards the clus- 
ters under study, we selected stars with spectral type earlier 
than AO. These have been identified from their location in 
apparent colour-colour (U — B versus B — V) and CM dia- 
grams, having V < 15 mag, and (B — V) < 0.70 mag, and for 
these stars, the MS spectral typ e have been obtained usin g 
UBV photometric Q-method l|johnson fc Morgan! Il953f ). 
The colour excesses are determined by subtracting intrin- 



sic colours from observed colours. The i ntrinsic colo u rs are 
derived from the colour relati on given bvlFitzGeraldl dl970h 
for (U - B) and (B — V); by I Johnson fc Morgan! (|l966f ) for 
(V - R) and (V - I) and by iKoornneej i| 19831 ) for (V - J), 
(V - H) and (V - K). To get the intrinsic colors for V - R 
and V — I, the present Cousins RI system a re conv e rted t o 
the Johnson RI using the relation given bv lBesselll l| 19791) . 
Table [5] lists the colour excess ratios estimated for both the 
clus ters alongwith the r atios for normal interstellar mat- 
ter l|Cardelli et all (|l989T l. The colour excess ratios generally 
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Table 4. Frequency distribution of the stars in the V, (V — I) diagram of the cluster and field regions. Nb, Ng and Nr denote the 
number of stars in a magnitude bin blueward. along and redward of the cluster sequence respectively. The number of stars in the field 
regions are corrected for area differences. Nc (difference between the Ng value of cluster and field regions) denotes the statistically 
expected number of cluster members in the corresponding magnitude bin. 
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Figure 6. The (U — B) versus (B — V) colour-colour diagram of 
the clusters. The continuous curve represents locus of Schmidt- 
Kaler's (1982) ZAMS for solar metallicity. The filled triangles are 
the emission type stars observed by us. 



Figure 7. The plot of (J — K) versus (V — K) colour-colour 
diagram of the cluster for the stars within the cluster radius. The 
solid line is the ZAMS of Z = 0.02 fitted for the marked values 
of the excesses while short dash lines show the error bars. 



agree within 3cr with those given for the normal interstellar 
extinction law. This indicates that the interstellar extinction 
law is normal in the direction of both the clusters. 

Using the relation R v = 1.1E(V - K)/E(B - V) 
ijWhittet fc van Bredalfl98Ch . we derived the values of Rv 
as 3.1 for NGC 637 and 3.2 for NGC 957 and this indicates 
normal interstellar extinction law in the direction of both 
the clusters. 



3.6 Distance to the clusters 

The ZAMS fitting procedure was used to derive distances 
to the clusters. Fig.[8]shows the intrinsic CMDs of probable 
cluster members as selected in < 13.2l for the clusters NGC 637 
and NGC 957. Employing mean E(B—V) for the clusters, we 
converted apparent magnitude and colours into the intrinsic 



one using the method as described in lYadav fc Sagarl (2002). 
A theoretical ZAMS (|Schaller et all 1 19921 ) for Z = 0.02 is 
visually fitted to the blue part of intrinsic CMDs in Vo, (U — 
B) ; Vo, (B- V) ; V , {V - R) and Vo, (V - I) . This 
gives an intrinsic distance modulus, (m — M)o, of 12.0 ± 
0.2 mag for NGC 637 and 11.7 ± 0.2 mag for NGC 957. 
The fact that we were able to find faint probable cluster 
members allow us to get a better definition of the cluster 
lower main sequence which in turn improves the estimation 
of the distances. The above distance moduli corresponds to 
a distance of 2.5 ± 0.2 kpc for NGC 637 and 2.2 ± 0.2 kpc 
for NGC 957. For NGC 637, it is similar to the value 2.5 



kpc derived by iHuestamendia et al.l (Il991 ) while it is les s 
than the value 2.75 kpc derived bv lPhelps fc Janesl (| 1994h . 
In case of NGC 957, our derived values are similar to the 
values given in literature (*JTJ). 
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Table 5. A comparison of extinction law in the direction of cluster with normal extinction law given bv lCardelli et al] l|l989h . 



Objects 


E(U~B) 
E(V-.J) 


E(B— V) 
E(V-J) 


E(V-R) 
E(V-J) 


E(V-I) 
E(V-.J) 


E{V-H) 
E(V-J) 


E{V-K) 
E(V-J) 


E(J-K) 
E(V~K) 


Normal value 


0.32 


0.43 


0.35 


0.70 


1.13 


1.21 


0.19 


NGC 637 


0.32±0.03 


0.45±0.04 


0.36±0.03 


0.69±0.05 


1.11 ±0.08 


1.19±0.07 


0.15±0.30 


NGC 957 


0.34±0.02 


0.47±0.01 


0.39±0.03 


0.75±0.05 


1.12 ±0.03 


1.19±0.02 


0.14±0.30 




-1 1 1 0.5 1 



(U-B)„ (B-V)„ (V-R) (V-I) 



Figur e 8. The intrinsic colour-magnitude diagram of the clusters. The continuous solid lines curve are the ZAMS given bv lSchaller et all 
(1992) for Z = 0.02 and the dotted lines are isochrones of log(age) = 6.8,7.0 and 7.2 (shown left to right respectively) for NGC 637 and 
NGC 957. 
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Figure 9. The K versus (J — K) and V versus (V — K) CMDs of the clu sters using probable cluster members. The solid curve represents 
the isochrones of log(age) = 7.0 for NGC 637 and NGC 957 taken from lSchaller et all dl992h for Z = 0.02. 



3.7 Age of the clusters 

Age of the cluster is determined by vis ually fitting the theo - 
retical stellar evolutionary isochrones (Sc haller et al.1 119921 ) 
for Z = 0.02 with the observed intrinsic CMDs of the clus- 
ters as shown in Fig. [5] We chose to use the isochrones of 
Z — 0.02, as the present target cluste rs are young and als o 
many recent studies, for example, bv lTwarog et al.1 (|l997h ■ 
does indicate that between 6.5 to 10 kpc galactocentric zone, 
the metallicity shows no gradient and its value is consistent 
with the s olar value within t he un certainties. Also for cluster 
NGC 637, iPhelps fc Janesl (|l994T ) considered solar metallic- 
ity. 

Isochrones of log(age) = 6.8, 7.0 and 7.2 are shown with 
dotted lines from left to right in Fig. [8] and these visually 
fitted isochrones to the brighter stars indicate that CMDs 
are represented well with an isochrone of about 10 Myr for 
both the clusters. Furthermore, the isochrones of log(age) 
= 6.8 and 7.2 Myr seem to confine blue and red end of the 
upper main sequence stars. The most massive stars in NGC 
637 and NGC 957 is about 13 Mq and it sets an upper limit 
of about 15 Myr for the MS age of the clusters. We have 
therefore adopted a mean age of 10 Myr with an uncertainty 
of about 5 Myr for both the clusters. 

The present age estimate for NGC 637 is not too 
different from the value s of 15 Myr as derived by 
IHuestamendia et all (Il99lf> and 4 Myr b y IPhelps fc Janesl 
( 1994 ). For NGC 957. lGerasimenkd (Il99ll ) has given an age 
of 5 Myr, which agrees within error with the present esti- 
mate. 

Using optical and near-IR data we re-determined dis- 
tance and age of both the clusters. We plot V versus (V — K) 
and K versus (J - K) CMDs in Fig [9] The theoretical 
isochrones given by ISchaller et ail (|l992l) for Z = 0.02 



and log(age) = 7.0 have been overplotted in the CMDs 
of NGC 637 and NGC 957. The apparent distance mod- 
uli (m — M) v ,(v-k) and (m — M)k,(j-k) turn out to be 
14.0±0.3 and 12.0±0.3 mag for the cluster NGC 637 and 
14.0±0.3 and 12.1±0.3 mag for the cluster NGC 957. Using 
the reddening values estimated in i)3.4l we derived a distance 
of 2.5±0.3 for NGC 637 and 2.2±0.3 kpc for NGC 957. Both 
age and distance determination for the clusters are thus in 
agreement with the estimates using optical data. However, 
scattering is larger due to the large errors in JHK mags. 

3.8 Gaps in the main-sequences 

A peculiarity that appears in the intrinsic CMDs is the pres- 
ence of gaps in the MS (see Fig. [SJ. A MS gap is present 
between 9 < V < 10 (11 < M/ Mp, < 8) in NGC 637 This 
gap has already been notic ed by IHuestamendia et al.1 (|l99ll ) 
and lPhelps fc Janesl l|l99St ). The deficiency of stars between 
11 < V < 12 (4 < M/Mq < 3) is also seen in NGC 957. 
To check the significance of the gap, we computed the prob- 
ability that a lack of sta rs in the mass interval is a re sult 
of random processes fsee lScald ll986l: iGiorgi et alJ liool) by 
using the relation P gap = (M sup /M inf ) ( -~ Nxx) , where M sup 
and Mi n f are the upper and lower masses of the gap and x is 
the exponent of mass function while N is the number of clus- 
ter members located above the gap. Adopting the MF slope 
in advance (gj), we find P 3ap = 4% for NGC 637 and 16% for 
NGC 957. The low resulting probabilities for both the clus- 
ters suggest that gaps in the MS are not due to the random 
processes. Hence, observed gaps in the MS of the clusters 
are real features, though we note that the significance test 
of gaps in MS are only approximate because the terms of 
the test were set after the data were observed. Specifically, 
upper and lower masses of the gap were derived after seeing 
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Figure 10. Slitless spectra of three emission line stars in NGC 
637 and NGC 957 



Figure 11. The normalised, continuum fitted spectra of the emis- 
sion star NGC 637-S1 in the wavelength range 3700-9000 A in the 
cluster NGC 637. 



the data. Thus, the quantitative statistical significance of 
the gaps has been overestimated. 

There are many reports about the MS gaps of other 
star clusters in t h e liter a ture jM crmilliod 1976; Forbes 1996; 
Herbst fe Milled 1 19821 ; iBaume et alj 120041 ; iKumar et all 
A detailed description 



20081 ). 



of 



;aps 



in several ope n 



cluste rs has recently been given by iRachford fc Canternal 
(2000). The explanation of the gaps found in open clus- 
ters ran ge from the onset of convection in the stellar en- 
velopes I Bohm-Vitense fc Canternal Il974h to peculiarities 
in the Balmer .j u mp and Balm er lines as suggested by 
iMermilliodl l|l976T ). |Ulrichlll971al lbh speculated that gaps in 
young clusters might be produced by the presence of 3 He iso- 
topes that halt gravitational contraction 1-2 mags above the 
main-sequence. Mermilliod's gaps occur in the range of B- 
type stars while Bohm-Vitense & Canterna's gaps are found 
at less massive stars. Observed gap in N GC 637 may be 
related to that found bv lMermilliodl (|l976Tl at the types Bl- 
B2. 



3.9 Emission type stars 

The Be star phenomenon in young open clusters is known 
to arise due to mass loss from early type stars during its 
evolution on the MS and are found to be maximum for clus- 
ters with 10-25 Myr age ijFabreeat fc Torrei6nll2000T ). About 
20% of the early-type (< B5) stars show Be phen omenon in 
young clusters of our Galaxy (Ma eder et al.lll999l ). However, 
in some young clusters, resulting from on-going/episodic 
star formation, an over-abundance of Be stars (> 30%) are 
caused due to presence of Herbig Be stars, which are usu- 
ally ass ociated with nebulosity and show large near-infrared 
excess (|Subramaniam et al.ll2005l . l2006f ). 

In the slit-less spectra, the emission-line stars are usu- 
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Figure 12. The normalised, continuum fitted spectra of the emis- 
sion stars in NGC 957 in the wavelength range 3700-9000 A. The 
top spectra corresponds to star NGC 957-S2 and the bottom one 
corresponds to star NGC 957-S1. 
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ally seen as continuum enhancement and they were identi- 
fied by visual inspection of each spectra. We could identify 
emission-line stars with a typical limit of about V ~ 16 mag 
in 10-min exposure. This roughly corresponds to a late F- 
type star for present clusters. The positions of the emission 
stars are identified from the _R-band images and the spectra 
are extracted from the slit-less spectral images. The spec- 
tra of individual stars are extracted with an aper t ure w idth 
of about 5 pixels (see ISubramaniam et al, 2005, 2006, for 
detailed description). We identified an emission line star in 
NGC 637 for the first time. We also identified two emission 
stars in NGC 957, which are previously known. The slitless 
spectra of one star in NGC 637 and two stars in NGC 957 
are shown in Fig. 10. These spectra cannot be used for quan- 
titative analysis, therefore, we obtained the slit spectra of 
the identified emission stars. The slit spectra of NGC 637-S1 
is shown in Fig. 1111 The spectra show H Q in emission, H/3 in 
emission within absorption. Two emission stars in NGC 957 
were also observed in the slit mode, the spectra are shown in 
Fig.[ll NGC 957-S1 is found to show H a , Hp and H 7 , along 
with a large number of Fell, OI and Paschen lines in emis- 
sion. NGC 957-S2 shows partially filled up Hp, faint emis- 
sions in Paschen lines along with H a emission. Photometric 
data with position of all three emission stars identified in 
the present study are listed in Table 6. Because of the CCD 
saturation in the present observations, photome tric data of 
two em ission stars in NGC 957 are taken from iHoag et al.l 
I 1961; and tabulated in Table 6. Thus, in the CMDs, the 
brighter and redder star is SI and the fainter and bluer star 
is S2. The values of H a equivalent width, stellar rotation, 
vsini, as estimated from the width of the Hel (447lA) line 
and the disk rotation, as estimated from the width of the 
H a profile are also tabulated. 

The spectral type of the star is estimated from its V 
and (B — V) magnitudes and using the estimated cluster 
reddening and distance. In NGC 637, the emission star is 
found to be a late B-type (B6.5V; NGC 637-S1) star with 
low infrared excess. This star is located well below the turn- 
off. The stellar vsini is found to be very large, 450 kms" 1 and 
disk vsini is found to be much smaller, 180 km s" 1 . Thus, this 
is a Classical Be star in the MS evolutionary phase and with 
large stellar rotation. Thus, this young clusters houses a fast 
rotating Classical Be star. 

The stars NGC 957-S1 and N GC 957-S2 in the cluster 
NGC 957 were earlier identified by iKohoutek fc Wehmeverl 
(1999) as emission type stars. The emission stars are found 
to belong to early B-type (SI - B2V and S2 - B2.5V) from 
photometry. These stars are located near the tip of the clus- 
ter MS, but below the turn-off. These stars have stellar vsini 
more than 300 kms" 1 indicating these are Classical Be can- 
didates. These stars also have low near-infrared excess which 
supports the above. The star with relatively large H a emis- 
sion, SI, is found to have a relatively slow rotating disk. 
Thus, this is another young cluster which houses Classical 
Be stars. 



4 LUMINOSITY AND MASS FUNCTION 

To construct the luminosity function for the cluster, we used 
the deepest V versus (V — /) CMD. Photometric criteria (see 
i|3.2p have been used to select the cluster members. The 
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Figur e 13. Mass function for NGC 637 and NGC 957 derived 
using Schallcr ct al. I Jl992f) isochrones. 



same envelope as made for selection of cluster members is 
also drawn in the V versus (V — I) CMD of the field region 
(see Fig[SJ. Stars are counted within this envelope, for both 
cluster and field region. Difference between the counts in 
two fields after accounting for the difference in area between 
the cluster and field regions will be the observed cluster 
luminosity function. To get the correct luminosity function 
it is very essential to derive completeness factor of CCD 
data. To know about the completeness of our photometric 
data, we performed experiments with artificial stars using 
ADDSTAR routine in DAOPH OT II. Detailed descrip tion 
about the experiment is given in lYadav fc Sagarl (|2002j) and 
ISagar fc Griffiths] (1998). The completeness factors derived 
in this way are tabulated in Table [7] for NGC 637 and NGC 
957 and it shows that the completeness factor is ~ 93% at 
20 mag for both clusters. The completeness factor of field 
region is assumed to be 100%. 

Using the cluster's paramete rs derived in th i s ana lysis 
and theoretical models given by ISchaller et all (|l992h we 
have converted LF to MF and the resulting MF is shown in 
Fig. 1131 The mass function slope can be derived by using 
the relation log|^j = —(1 + a;) x log(M)+constant, where 
dN represents the number of stars in a mass bin dM with 
central mass M and x is the slope of MF. The slope of 
the MF is x = 1.65 ± 0.20 and 1.31 ± 0.50 for the cluster 
NGC 637 and NGC 957 respectively. Our derived values of 
mass function slope is in agreement wit hin the error with 
the value 1.35 given bv ISalpeterl (|l955h for field stars in 
solar neighbourhood. 

There are many MF studies available in the litera- 
ture using open star clusters in the Milky Way. Recently, 
IPiskunov et al.l (|2004h studied 5 young open star clusters 
and found that stellar mass spectra of these clusters are 
well represented with a power law very similar to Salpeter 
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Table 6. Photometric information about emission stars. Photometric data of NGC 957-S1 and NGC 957-S2 star are taken from Hoag 
et al. (1961). 



Objects 


Coordinate (J2000) 


V 


(U-B) 


(B-V) 


( V - R) 


(V-J) 


H Q EW 


vsini(disk) 


vsini(star) 


NGC 637-S1 


a = 0l'M3 r ™23 s .l 


13.67 


-0.03 


0.55 


0.38 


0.78 


-21.97 


179.65 


450.37 




S = 64°0l'l9".2 


















NGC 957-S1 


a = 02 fe 33" I 08 s .2 


11.13 


-0.38 


0.66 






-37.99 


184.32 


342.72 




S = 57° 28' 10". 8 


















NGC 957-S2 


a = 02 h 33"*10 s .5 


11.99 


-0.28 


0.54 






-16.96 


251.62 


308.85 




S = 57°32'53".9 



















Table 7. Variation of completeness factor (CF) in the V, (V — I) 
diagram with the MS brightness. 



V mag range NGC 637 NGC 957 
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value within the uncertainties. MF study of iPhelps fe Janes! 
l| 19931 ) using a sample of seven star clusters also conclude 
the Salpeter type MF slope in these clusters. These studies 
alongwith the others ([Sanner fc Geffertl l200ll ; ISagar et alj 
l200li ; lYadav fc Sagarll2002l . |2004| ) found their results consis- 
tent with the Salpeter value. 



5 MASS SEGREGATION 

In order to study the mass segregation effect in the clusters, 
we plot the cumulative radial stellar distribution of stars for 
different masses in Fig [14] and mass segregation effect is seen 
for both the clusters, meaning, higher-mass stars gradually 
sink towards the cluster center than the lower mass stars. 
Further, We performed K-S test to see the statistical sig- 
nificance of mass segregation. For NGC 637, we divided the 
MS stars in three magnitude range i.e. 12.0 < V < 15.0, 
15.0 < V < 18.0 and 18.0 < V < 20.0 corresponding to the 
mass range of 9.0 < M/M e < 2.5, 2.5 < M/M Q < 1.2 and 
1.2< M/M Q <0.8. For NGC 957, the magnitude range are 
12.0 < V < 14.0, 14.0 < V < 16.0 and 16.0 < V < 18.0 
corresponding to the mass range of 8.0 < M/Mq < 3.5, 
3.5 < M/Mq < 2.0 and 2.0 < M/M e < 1.0. The K-S test 
shows the mass segregation effect at confidence level of 99% 
for NGC 637 and 70% for NGC 957. 

Mass segregation effect can be due to dynamical evolu- 
tion or imprint of star formation or both. In the lifetime of 
star clusters, encounters between its member stars gradually 
lead to an increased degree of energy equipartition through- 
out the clusters. In this process the higher mass cluster mem- 
bers gradually sink towards the cluster center and transfer 
their kinetic energy to the more numerous lower-mass stel- 
lar component, thus leading to mass segregation. The time 
scale on which a cluster will have lost all traces of its ini- 
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Figure 14. The cumulative radial distribution of stars in various 
mass range. 



tial conditions is well represented by its relaxation time Te, 
which is given by 

8.9 x W 5 N 1/2 R 3 h /2 
E ~ < m >V2 log(0AN) ( ' 

where is the number of cluster members, Rh is the half- 
mass radius of the cluster and < m > is the mean mass 
of the cluster stars (Spi tzer fc Hartlll97ll ). The value of Ru 
has been assumed as half of the cluster radius derived by 
us. Using the above relation we estimated the dynamical 
relaxation time T E = 10 Myr for NGC 637 and NGC 957. 
This indicates that the age of both clusters is same as its 
relaxation age. Therefore, we conclude that both the clusters 
are dynamically relaxed. 



6 CONCLUSIONS 

We studied two open star clusters NGC 637 and NGC 957 
using UBVRI CCD and 2MASS JHK data. Results ob- 
tained in the analysis are the following. 

(i) The radii of the clusters are obtained as 4!2 and 4.'3 
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which corresponds to 3.0 and 2.8 pc respectively at the dis- 
tance of the clusters NGC 637 and NGC 957. 

(ii) From the two colour diagram, we estimated E(B — 
V) = 0.64 ± 0.05 mag for NGC 637 and 0.71±0.05 mag for 
NGC 957. The JHK data in combination with the optical 
data provide E(J - K) = 0.23±0.20 mag and 0.23±0.20 
mag while E(V - K) = 1.70±0.20 mag and E(V — K) = 
2.00±0.20 mag for NGC 637 and NGC 957 respectively. 
Analysis indicates that interstellar extinction law is normal 
towards the direction of both the clusters. 

(iii) Distances to the cluster NGC 637 and NGC 957 are 
determined as 2.5±0.2 and 2.2±0.2 kpc respectively. These 
distances are supported by the distance values derived us- 
ing optical and near-IR data. An age of 10±5 Myr is deter- 
mined for both NGC 637 and NGC 957 by compari ng the 
isochrones of Z = 0.02 given bv lSchaller et all (|l992T l. 

(iv) We discovered one emission type star in NGC 637. 
Two emission type stars in NGC 957 are also confirmed by 
the present study. These are classical Be star candidates. 
The two emission stars in NGC 957 are about to leave the 
MS while the emission star in NGC 637 is in the MS evolu- 
tionary phase. 

(v) The mass function slopes x — 1.65 ± 0.20 and 1.31 ± 
0.50 are derived for NGC 637 and NGC 957 by considering 
the corrections of field star contamination and data incom- 
pleteness. Our analysis indicate that both the cluster are 
dynamically relaxed and one plausible reason of this relax- 
ation may be the dynamical evolution of the clusters. 
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